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Abstract. The model by Braginskii for the viscous stress tensor is used to 
determine the shear and gyroviscous forces acting within a toroidally confined plasma. 
Comparison is made to a previous evaluation which contains an inconsistent treatment 
of the radial derivative and neglects the effect of the pitch angle. Parallel viscosity 
contributes a radial shear viscous force which may develop for sufficient vertical 
asymmetry to the ion velocity profile. An evaluation is performed of this radial viscous 
force for a tokamak near equilibrium which indicates qualitative agreement between 
theory and measurement for impure plasma discharges with strong toroidal flow. 
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1. Introduction 

The evaluation of tlie Braginskii model of plasma viscosity pLj for toroidal confinement is 
hampered by the complicated geometry. Previous attempts [3l H] have not treated the 
radial derivative in a consistent manner, nor have they fully taken into account the effect 
of the magnetic field pitch angle on the viscous stress tensor. As the determination of 
this quantity is of intrinsic interest to the plasma theorist as well as to fusion engineers, 
here we reevaluate the shear and gyroviscous forces for a toroidal plasma. 

This paper focuses on the viscosity model derived by Braginskii for a large flow 
velocity V ~ Vth and assumes a ratio of gyrofrequency to self-collision frequency 
uj/u = LOT ^ 1. We note the recent derivations by Mikhailovskii and Tsypin [5j, 
Catto and Simakov [HITIIH], and Ramos [9] consider an ordering with small flow velocity 
V <^ Vth which introduces terms dependent on the heat flux. We will neglect such effects 
as well as any arising from trapped particles [H [TOl |Tl], [12] as a first approximation to the 
extension of the classical Braginskii viscosity to toroidal geometry. In situations where 
the assumptions of applicability are not met, these additional terms may be necessary. 

After working through in detail a prevailing derivation of the toroidal gyroviscous 
torque, we consider an alternative approach based on the geometric nature of 
tensor quantities. For both models we will use the concentric, circular flux surface 
approximation; however, the treatments of the radial derivative will be distinct. 
Restricting consideration to flux surfaces with constant density yet allowing for vertical 
asymmetry to the ion toroidal velocity profile yields a radial shear viscous force which 
should be accounted for when interpreting experimental measurements [131 IE]- An 
evaluation from data collected at DIII-D [l5j indicates qualitative agreement between 
theory and measurement for impure plasma discharges with significant toroidal rotation. 

2. Stacey-Sigmar model for gyroviscous torque 

We begin by considering in detail the derivation of the gyroviscous torque as given 
by Stacey and Sigmar [3J. Using a concentric, circular toroidal flux surface geometry 
{r,6,(j)) with major radius Rq allows one to write R = Rq + rcos6 = Ro{l + ecosO) 
and B = B°/(l + ecos6') with V ■ B = to set 5^ = 0. Derivatives d/dx may 
be written dx, and we assume toroidal symmetry 0. Considering a single ion 

species, the density on the flux surface is expanded about a radially dependent constant 
as n = n°(r)[l + e{n'^cos6 + n*sin^)], with similar expressions for the components 
of the fluid velocity V, and the radial dependence dr = Ro^d^ of the poloidal 
variations are neglected when taking the radial derivative of an expanded quantity, 
drU ~ {n/n^)drn^ , despite the exphcit appearance of the factor e as well as any radial 
dependence of the poloidal coefficients n'^'**(r). The poloidal derivative is taken explicitly, 
den = n^e{n^ cos 6 — n'^sinO). In this section, the thermal energy ksT — > T = 
assumes no poloidal variation so that the pressure takes on the dependence of the density, 
p = nT = p^{n/n^), for consistency of comparison. One assumes here that both the 
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ions and the electrons are isothermal. 

For an equilibrium density dtU — > the continuity equation reads h = dn/dt = 
V ■ nV = {drvRnVr + deRnVg)/rR, where n is the particle source rate. The unity 



moment flux surface average (A)^ = de{l + e cos 9) A{e) /27r yields the differential 
equation n° = n^V^/r + nPdrV^ + V^driiP which has nontrivial solution for the radial 
velocity only with nonzero source rate; for n 0, Vr — > 0. Then the cosine 
(A)^ = ^^'^ decos9{l + ecos9)A{e)/2n and sine {A)g moments relate the density 

— c — s 

variations to those of the poloidal velocity, Vo — —rf — 1 and Vq — —n^. The flux 
surface averaged gyroviscous torque is derived from the rate of shear tensor under the 
assumption that \B^/B\ 1 and \Bq/B\ fa as 

nTm/zeB is the gyroviscous coefficient. We spht the integrand into five 
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collecting their units into the common factor lu = nPTmV^ /2zeB^Ro — rflOP /2, and 
evaluate their leading contributions thusly, 
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where L]^^ = — 9,. lnX° is the inverse radial gradient scale length for X, giving a net 
gyroviscous torque of 

R^-V-U^'^-"^''^' 
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{r (J2 Lx') [V (4 + n'^) +n^{l- n')] - 4T^/ -n'} , (12) 



where ^L^^ = + Ly^ — . The Stacey-Sigmar model is seen to correspond 
to the contribution from 7i + 74. The evaluation of I2 and arises from considering 
the poloidal dependence of it! and B as geometrically induced rather than treating 
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them as expanded quantities, and the term I3 represents the change in the gyroviscous 
coefficient dependent upon the gyrofrequency, as drB^ = {B^ / B^)drBQ when 9,. 5° = 0. 
Interestingly, the final terms in ( |T2l) are of opposite sign, hence would tend to cancel, 
those of the Stacey-Sigmar model. The derivation for the poloidal shear viscosity [21 H] 
proceeds similarly with the same caveats and pitfalls. 



3. Viscous force for toroidal confinement 



Next we consider an alternate derivation making use of the transformation properties 
of the viscosity tensor. We also will treat the radial derivative consistently throughout. 



While our expansion looks the same. 
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specify the poloidal coefficients to have no radial dependence, dj.n'^''^ = 0, and the radial 
derivative is explicitly drti = {n/n^)drn^ + n^{n/'nP — l)/r. The unity moment of the 
continuity equation remains the same, and remembering that physically it is the dynamic 
equation for the density, we use its cosine and sine moments to write n'^ = —Vq — 1 
and n'^ = —Ve when Vr That we have used its unity moment to determine V^: 

needs to be recalled when addressing the radial component of the equation of motion. 
In the presence of a radial fiow, the density coefficients acquire a correction factor 
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Y(l + 2r9rln?2°) which would introduce an r dependence unless 
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The poloidal dependence of the ion temperature is addressed through the equilibrium 
equation of state V ■ Vp + 7pV ■ V = for adiabatic index 7, with solution to 0(e) of 
T^'' = %nf' for 7i = 7i - 1. 

Moving on, the rate of shear tensor in general is written as 
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and reduces for toroidally symmetric (r, 6, 
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which is a geometric object independent of one's choice of coordinate system. Identifying 
(r, 6, (f)) as the primed coordinate system, the transformation to an unprimed coordinate 
system {x,y,z) is given by W = G^WG, where G G S0(3) is a general parity 
preserving rotation matrix in M.^. Braginskii's theory is derived for coordinate axes 
orientated to the magnetic fleld, thus for {x,y,z) = (r, ±, ||) where \\ = b = B^/B^ = 
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as given in Figure [H Our strategy is to take W W — 11 — >■ 11', where 11' = GIIG"^. 
The viscosity coefficients are in the ordering t]q ^ t]4 ^ ri2 for a strong magnetic field 
ujT ^ 1, and we will neglect the perpendicular contribution 771,772 0. The shear 
viscosity reduces to = — "^oWq, and the gyroviscosity to YIq = "^4(^3/2 + W4), 
where 770 — nTr and 774 ~ nTm/zeB such that f]o/Vo = (T/T°)^/^ and rjA^/rfl = 
{n/n^){T/T^){R/Ro). Expficitly, 
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In our evaluation of the viscosity, all derivatives are taken exphcitly, neglecting only 
the poloidal dependence of the Coulomb logarithm. The components of the viscous force 
are given by (V ■ 11)^ = Y^/sl^is^l^a + Y.'y(X^'y^i3'y + ^p-y^'ya)], and we write the toroidal 
viscous torque 
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noting there is no 770 contribution to Ilr<t>- Evaluating its leading components, we find 
that the shear viscous torque vanishes, {Rep ■ V ■ II'^.) ~ 0, and the gyroviscous torque is 
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whereupon restricting consideration to fiux surfaces with constant density, = and 
Ve = —1, denoted by the symbol yields a non-vanishing torque {R(f) ■ V ■ H'q) =^ 
-A{riyRo)V%M + '^blbj - bp^)^' ■ The poloidal viscous force is 
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where Uro has only an 774 component, and its shear contribution evaluates to 
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The radial viscous force is 
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where we have used drB^ = ^iqJ^ — Bg/r. According to this theory, the poloidal velocity 
must go as a quadratic (or greater) near the magnetic axis, ~ as r — > 0, else the 
radial gyroviscous force will diverge. An order of magnitude estimate indicates that this 
force is on par with the radial inertial force, (nm(V ■ V)V) ^ — n°mV^g°^/r ~ lONt/m^. 
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4. Radial force from vertical asymmetry 

4-1. Motivation 

Returning to the radial shear viscous force for constant density flux surfaces (|38|) . we 
see that a radial viscous force may develop from an interplay between the magnetic 

s 

field pitch angle 6e 7^ and a vertical asymmetry to the toroidal velocity 7^ 0. 
We are not aware of this result appearing previously in the literature, as it disappears 
for an ordering hg -C h^. Considering now the effects of a single null vertical divertor 
configuration as commonly found in a tokamak discharge, we expect there will exist a 
vertical asymmetry to the plasma, particularly to the toroidal velocity. Thus, the radial 
shear viscous force should make a contribution to the radial equation of motion, 

- (f ■ [nm (V ■ V) V + V • Us]) - {drP - nze {B^Ve - BeV^)) (46) 
n^mVf/r -{r-V -Us) ~ n^zeE^'^P , (47) 

where the right hand side is commonly identified as the radial electrostatic field yet 
consists actually of experimental measurements of density, temperature, magnetic field, 

s 

and impurity velocity. Thus, we are led to a model with one free parameter, , 
which we can compare to data taken from a tokamak. As has been determined 
from the equilibrium continuity equation, the unity moment of the radial equation 
of motion becomes the appropriate place to make a comparison between theory and 
experiment. We mention here that, for the single-fluid model, the natural degrees of 
freedom associated with the conservation equations of mass, energy, and momentum are 
?T,(V), T{n,\'), and Y{n,T), respectively. 

' — ^ 

Restricting the range of the parameter to the region \emaxV(t> \ = \ V^\ < 1 keeps the 
expansion for the toroidal velocity from changing sign, where Emax = a/Ro for minor 
radius a. For a geometry with 6 downward at the outer midplane, a positive/negative 
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Table 1. Best fit parameters. 



shot 




a 


Fig. [7] 


119306 


0.29 


15.77 


(a) 


119307 


0.50 


14.11 


(d) 


121433 


0.83 


16.44 


(b) 


121453 


1.00 


84.50 


(e) 


122338 


1.00 


49.51 


(c) 



value of indicates an increase/decrease in the toroidal velocity at the poloidal null, 
with an opposite effect at the upper vertical midplane. Minimizing x^, the sum of 
squared residuals, over yields the parameter of best fit as well as the normalized 
root-mean-square residual a. 

4-2. Data Analysis 

To evaluate the comparison of the predicted radial viscous force with experimental 
measurements, we need the profiles for density, temperature, poloidal and toroidal 
magnetic field, and the rotation profiles for at least one species of ion. To analyze 
discharges at DIII-D [15], we retrieve the necessary magnetic geometry and pressure 
profiles from EFIT [16], and the rotation profiles for Cg are corrected for deuterium's 
energy dependent charge exchange cross section [13j. Collision rates are given by the 
NRL Formulary [T7|. The measurement profiles for our selection of shots are given in 
Figures [2] through [6l Times given are representative of the available measurement times, 
and the shots 121433 and 121453 have had helium and neon impurities injected. 

Our results presented in Figure [7| compare the prediction on the left hand side of 
( H6l) with the measurements on the right hand side, where the Theory values include the 
inertial term. For three of the five shots considered, the theoretical force density can 
agree with the measurement for a fitted parameter | V"^ | < 1 , found in Table H], where the 
positive values of the parameter indicate that the plasma increases in toroidal velocity 
near the poloidal null and a is normalized to the units of the figure. However, for the 
other two, the maximum of likelihood lies outside the prior range and the posterior 
is correspondingly truncated. (A model with all 4 velocity parameters flM|) -(37). thus 
including the density and temperature poloidal variations as related above, did not 
noticeably improve the quality of fit, hence is discarded by Ockham's principle |18j.) 
For shot 121453, the theory cannot reach the magnitude of the measurement profile on 
account of the low values of the toroidal velocity and ion temperature, though its shape 
is consistent with the others, and for shot 122338 the shape of the profile cannot be 
matched by the present theory. Were it not for these two shots, we would conclude 
that the Braginskii theory of viscosity is entirely consistent with the experimental 
measurements. 
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5. Discussion and Conclusion 

How much impact has our neglect of any heat flux had on our results? The heat 
flux contribution to the generalized viscosity appears at orders consistent with the 
perpendicular and gyroviscosity [1] for a strong flow, though recent authors suggest 
it contributes to the parallel viscosity as well for subsonic flow. For a model based on 
an equilibrium evaluation, however, the temperature is assumed to be steady in time, 
hence no heat flows on the timescales of interest. Such assumption is consistent with 
our assumption of constant density hence constant temperature on the flux surfaces. 
The parameter is independent of Vq, thus also of and T'^. Our results suggest 
that a viscous force may develop in conjunction with vertical asymmetry introduced by 
a divertor conflguration which has not been considered previously in the literature. 

Whenever one compares a theory to a measurement, one must consider whether 
the experiment has actualized the assumptions of applicability. In a real tokamak, a 
multitude of processes may be occuring which are not considered here, and the failure 
of the model to account for the radial force density proflle of two of the flve shots may 
reflect these neglected phenomena. The evaluation of the Braginskii viscous force is 
straightforward using axes rotation, and the introduction of higher order effects does 
not obviously detract from the parallel contribution. Certainly improvements to the 
geometry to handle flux surfaces of arbitrary cross sectional shape can be considered [19]; 
however, the major discrepancy of the concentric circular flux surface approximation is 
in the horizontal direction, not the vertical direction driving the radial viscous force. 
Nevertheless, the qualitative agreement found for the majority of the shots considered 
here suggests the radial shear viscous force is an effect worthy of further investigation. 

The Braginskii model for plasma viscosity in a strong magnetic fleld has been 
evaluated in the concentric circular flux surface approximation for toroidal confinement 
retaining the parallel and gyroviscous contributions by means of the rotation matrix 
between coordinate axes in the lab frame and those aligned to the magnetic field. 
Explicit attention is paid to the radial derivative of expanded quantities. Restriction to 
fiux surfaces of constant density yields a model allowing for vertical asymmetry in the 
impurity ion toroidal velocity which drives a radial shear viscous force with a magnitude 
on the order of that attributed to the radial electric field. A prediction for the vertical 
asymmetry to the toroidal velocity obtains from fitting the radial shear viscous force 
to the experimentally measured radial force density. An evaluation from data indicates 
qualitative agreement between the theoretical and measured force densities for several 
discharges with strong toroidal fiow. 
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Figure 1. (Color online.) The transformation of vectors and tensors between 
orthonormal coordinate axes is given by the rotation matrix, which may be found 
from the components of one set of unit vectors along the other set of coordinate axes. 
Here, y = b^O — hgcj) and z = bgO + b^cf). 




normalized radius r/a 



Figure 2. (Color online.) Measured profiles for shot 119306 at time 2955ms. 
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normalized radius 



Figure 3. (Color online.) Measured profiles for shot 119307 at time 2205ms. 
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Figure 4. (Color online.) Measured profiles for shot 121433 at time 3075ms. He and 
Ne also were present. 



14 




Figure 5. (Color online.) Measured profiles for shot 121453 at time 3075nis. He 
and Ne also were present. Note the low value for the toroidal velocity (c) and ion 
temperature (f). 
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normalized radius 



Figure 6. (Color online.) Measured profiles for shot 122338 at time 2750ms. 




Figure 7. (Color online.) Comparison of theory (x) to measurement (0) in (HH). The 
parameters of best fit are in Table [Tl 



